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Abstract 

Following the formation of an excited muonic atom, inner shell transitions 
may proceed without photon emission by inverse internal conversion, i.e. the 
muonic excitation energy is transferred to the nucleus. In actinides, the 2p — > Is 
and the 3d — > Is muonic transitions result in excitation of the nuclear giant 
dipole and giant quadrupole resonances, respectively, which act as doorway 
states for fission. The nuclear excitation energy is typically 6.5 — 10 MeV. Be- 
cause the muon lifetime is long compared to the timescale of prompt nuclear 
fission, the motion of the muon in the Coulomb field of the fissioning nucleus 
may be utilized to learn about the dynamics of fission. If there is large friction 
between the outer fission barrier and the scission point the muon will remain 
in the lowest molecular energy level and emerge in the Is bound state of the 
heavy fission fragment. On the other hand, if friction is small (i.e. the nuclear 
collective motion is fast) there is a nonvanishing probability that the muon may 
be promoted to higher-lying molecular orbitals, e.g. the 2pa level, from where 
it will end up attached to the light fission fragment. Therefore, theoretical 
studies of the muon- attachment probability to the light fission fragment, Pl, 
in combination with experimental data can be utilized to analyze the dynamics 
of fission, and nuclear energy dissipation in particular. In this way, one may be 
able to distinguish between nuclear energy dissipation arising from two-body 
collisions and "one-body friction" caused by the mean field. We study the dy- 
namics of a muon bound to a fissioning actinide nucleus in two different ways: 
first, we solve the non-relativistic time-dependent Schrodinger equation using 
the Born-Oppenheimer expansion method in terms of quasimolccular wavcfunc- 
tions. In the second approach, we solve the relativistic time-dependent Dirac 
equation on a 3-D Cartesian lattice utilizing the B-Spline collocation method. 
For 2 |jNp we find a dissipated energy of order 0—10 MeV and a fission time 
delay due to friction of up to 2 • 10~ 21 s. 



1 



Muon-induced fission 



2 



1 Introduction 

Nuclear physics experiments with \i~ beams provide information on fundamen- 
tal symmetries and interactions. Muonic atoms, in particular, have proven ex- 
tremely useful in examining the electromagnetic properties of nuclei, e.g. electric 
charge distributions and multipole moments, because the muon has a high posi- 
tion probability density inside the nucleus owing to its small Compton wavelength 
A c = h/(mu,c) = 1.87 fm pQ. After muons have been captured into high-lying single 
particle states they form an excited muonic atom. Inner shell transitions may pro- 
ceed without photon emission by inverse internal conversion |2j, i.e. the muonic exci- 
tation energy is transferred to the nucleus. In actinides, the El : (2p — > Is, 6.6MeV) 
and the E2 : (3d — > ls,9.5MeV) muonic transitions result in excitation of the nu- 
clear giant dipole and giant quadrupole resonances, respectively, which act as door- 
way states for fission. Following these atomic transitions, the muon will be in the 
ground state of the actinide atom and can be utilized to probe the fission dynamics. 
Finally though, the muon is going to be captured by one of the fission fragments. 
However, nuclear capture as a result of the weak interaction occurs on a time scale 
of order 10 -7 s which is many orders of magnitude larger than the time scale of 
fission. 

From a theoretical point of view, prompt muon-induced fission has several at- 
tractive features. Because the nuclear excitation energy equals or exceeds the fission 
barrier height it is permissible to treat the fission dynamics classically (no barrier 
tunneling). The muon dynamics is determined by the electromagnetic interaction 
which is precisely known; hence, the process can be calculated, at least in principle, 
with any desired precision. Our main task is the solution of the Dirac equation 
for the muon in the presence of a time-dependent external Coulomb field which is 
generated by the fission fragments in motion. 

We will demonstrate that the muon attachment to the light fission fragment 
depends on the nuclear friction between the outer fission barrier and the scission 
point. In this context, nuclear friction is defined as the irreversible flow of en- 
ergy (and linear or angular momentum) from collective to intrinsic single-particle 
motion We include in our classical dynamical calculations for the fission mode a 
linear friction force to account for energy dissipation via neutron and photon emis- 
sion. Through muon-induced fission one expects to gain a deeper understanding of 
the energy dissipation mechanism in large- amplitude nuclear collective motion. 

The prompt muon-induced fission process is most easily understood via a "cor- 
relation diagram" , i.e. one plots the single-particle energies of the transient muonic 
molecule as a function of the internuclear distance (see Fig. 3 of ref.[l]). If there is 
a large amount of friction during the motion from the outer fission barrier to the 
scission point the muon will remain in the lowest molecular energy level lsa and 
emerge in the Is bound state of the heavy fission fragment. If, on the other hand, 
friction is small and hence the nuclear collective motion is relatively fast there is a 
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nonvanishing probability that the muon may be promoted to higher-lying molecular 
orbitals, e.g. the 2pa level, from where it will end up attached to the light fission 
fragment. Therefore, theoretical studies of the muon-attachment probability to the 
light fission fragment, Pl, in combination with experimental data can be utilized to 
analyze the dynamics of fission, and nuclear energy dissipation in particular. 

There are two different mechanisms that contribute to nuclear energy dissipa- 
tion: two-body collisions and "one-body friction". The latter is caused by the 
moving walls of the self-consistent nuclear mean field. The role played by these 
two dissipation mechanisms in fission and heavy-ion reactions is not yet completely 
understood. In 1976, in a pioneering work Davies, Sierk and Nix [5] calculated the 
effect of viscosity on the dynamics of fission. Assuming that friction is caused by 
two-body collisions they extracted a viscosity coefficient /x = 0.015 Tera Poise from 
a comparison of theoretical and experimental values for the kinetic energies of fission 
fragments. The corresponding time delay for the nuclear motion from the saddle to 
the scission point was found to be of order At = 1 x 10 -21 s. However, in one-body 
dissipation models the time delay is an order of magnitude larger. 

Theoretical studies in combination with rather scarce experimental data on 
prompt fission obtained so far give rise to a definite answer: Estimates of the 
energy dissipated in the case of the one-body mechanism were made in the two- 
center harmonic oscillator model in analytical form j^j. For this purpose, the Born- 
Oppenheimer expansion method developed for calculation of the muon promotion 
probability, was applied for assessment of the nucleon promotion rate to higher 
quasimolecular orbitals, arising from the translation of the nascing fragments. The 
last process provides the microscopic mechanism for the one-body dissipation. The 
obtained estimate of the dissipated energy, Edi ss < lMeV, allows one to rule out the 
one-body mechanism, even with all the appoximations of the model kept in mind. 

Several experimental techniques are sensitive to the energy dissipation in nuclear 
fission. At high excitation energy, the multiplicity of pre-scission neutrons [7] or 
photons depends on the dissipation strength. At low excitation energy, the 
process of prompt muon-induced fission |U] provides a suitable "clock" . This process 
will be discussed here. 

Nonrelativistic calculations for muon-induced fission have been carried out by 
several theory groups, starting in 1980 [13 E3 ED E] j but relativistic calculations 
using the Dirac equation for the muon only became feasible in 1992 |14l lH [To1ll6j . We 
discuss here in detail our theoretical approach and the numerical implementation, 
and we compare our results with experimental data obtained at the Los Alamos 
Meson Physics Facility (LAMPF) [T71 Q2], at the Tri-University Meson Facility 
(TRIUMF) H2| ED|, and at CERN and the Paul Scherrer Institute (PSI) [2UE21I2S1 
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2 Prompt and Delayed Fission Induced by Muons 

Following the irradiation of a target with a [x~ beam the muons lose most of their 
kinetic energy by ionization in the target material within 10 -9 to 10 _10 s. Once 
their velocity has become comparable to the orbital electron velocities characteristic 
of these atoms, they are slowed down further by inelastic collisions with valence 
electrons and are finally captured into high-lying states (Ua ~ 14) forming a muonic 
atom. The theoretical aspects of the interaction of muons with condensed matter 
were first studied by Fermi and Teller j^Hj and were later explored in more detail by 
Wu and Wilets [22] an d by Kim . Because all muonic bound states are unoccupied 
the muon will cascade down to the ground state within 10~ 13 s. From the outer shells 
the excited muonic atom decays preferentially by emission of Auger electrons. Since 
AE increases rapidly for the inner shells, the transitions between levels with n < 5 
are dominated by mesic X-rays. Alternatively, the transitions may proceed without 
emission of radiation via inverse internal conversion. From the K-shell, the muon 
disappears at a characteristic rate A = Ao + A c , where Ao = (2.2 x 10 _6 s)~ 1 denotes 
the free leptonic decay rate and A c the nuclear capture rate; A c depends upon the 
charge and mass of the nucleus (Goulard-Primakoff formula 28 ) and is of order 
(7.5 x 10 _8 s) _1 for actinides JH]- Muons stopped in an actinide target may induce 
nuclear fission in two different ways: 

• Delayed Fission Following Nuclear Muon Capture 

The muon is captured by a proton inside the nucleus and forms a neutron and 
a muon neutrino 

f i- + (Z,A)->(Z-l,A)* + v ll . (1) 

Even though most of the energy is taken away by the neutrino, the average 
nuclear excitation energy is 15 — 20MeV which is well above the fission barrier 
for actinides Et = 5 — 6MeV . Fission via nuclear muon capture is delayed, i.e. 
it occurs with the characteristic mean lifetime of the weak decay process, r cap t = 
(7-8)xl(T 8 s. 

• Prompt Fission Resulting From Inverse Internal Conversion In Muonic Atoms 
In this case, the excitation energy of the muonic atom is transferred to the 

nucleus by an internal conversion process (nonradiative transition) and the muon 
ends up in the K-shell of the muonic atom 

(Z,A)(/i-)*->(Z,A)V-. (2) 

For the innermost atomic transitions in an actinide muonic atom, the transition 
energy generally exceeds the fission barrier height. The result is prompt fission 
in the presence of the muon, since the muon is not annihilated by this process, 
in contrast to fission resulting from nuclear muon capture. The nucleus will be 
surrounded by the muon during the entire fission process, unless the muon is ionized. 
Eventually, the muon will decay by nuclear muon capture from the fission fragments. 
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Experimentally, both fission modes can be distinguished because of their different 
time scale. In this paper, we focus on prompt muon-induced fission. This process was 
first discussed by Wheeler |29j and considered in more detail by Zaretski and Novikov 
|3U| . It is important to know the specific atomic transitions that are responsible for 
prompt fission. Fig.l of ref.jl] shows the Coulomb interaction energy between the 
muon and a 238 U nucleus as well as the binding energies of the lowest bound states. 
Even though EO transitions such as 2s — ► Is and 3s — ► Is exhibit the largest internal 
conversion rates they do not contribute to fission because they lead to excitation 
of the giant monopole resonance which is spherically symmetric and much too high 
in energy. On the other hand, the {El : 2p — > Is) and the {E2 : 3d — > Is) 
transitions result in excitation of the electric giant dipole and quadrupole resonances, 
respectively, both of which act as doorway states for fission. Let us consider the 
specific case of 238 U: The giant dipole resonance is located at Egdr = 12.8MeV 
and has a width T = 6MeV^£Q; for the T = giant quadrupole resonance the 
corresponding numbers are Eqqr = 9.9MeV and T = 6.8MeV^^. According to 
Teller and Weiss [33\ and Karpeshin and Nesterenko [2] it is very probable that the 
3d — > Is radiationless transition will be dominant for muon-induced fission, because 
its transition energy of 9.6MeV is very close to the peak of the giant quadrupole 
resonance whereas the 2p — * Is transition energy of 6.6MeV is far off the center of 
the giant dipole resonance. Experimentally, the situation is controversial: Johansson 
et al. |22| measured muonic X-rays in coincidence with prompt fission in 238 U. From 
the muonic X-ray intensity ratios for prompt and delayed fission they conclude that 
(74 ± 15)% of all prompt events can be attributed to the 3d — > Is radiationless 
transition and only (26 ± 15)% to the 2p — * Is transition. On the other hand, 
Kaplan et al. j2Hj find in similar studies a predominance of the El transitions in 
their prompt fission data. 

As mentioned in the Introduction, if the nuclear motion is relatively fast (low 
friction) there is a nonvanishing probability that the muon may be promoted to 
higher-lying molecular orbitals, e.g. the 2pa level from where it may end up at- 
tached to the light fission fragment. This gave rise to hopes that the muon attach- 
ment probability to the light fragment Pl would depend on the pre-fission scenario. 
In this case, theoretical studies of the muon-attachment probability to the light fis- 
sion fragment in combination with experimental data can be utilized to probe the 
dynamics of prompt fission. 

On the other hand, in a series of papers it was pointed out that the fragment 
motion occurs quasiclassically. The fragment velocity exceeds by ~ 6 times that of 
the muon in the orbit. Moreover, it satisfies the Massey criterium of adiabaticity 
which does not promise any significant transition probability. 

Fortunately, the presence of the avoided crossing of the lscr and 2pa levels favors 
the transition. As is known, in this case the transition occurs in a narrow vicinity of 
the crossing. At the same time, the pseudocrossing occurs well beyond the scission 
point, at a distance which is twice as large as that of the neck rupture for the most 
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likely pairs of fragments. At this distance, the trajectory of the relative fragment 
motion R(t) is governed by their mutual Coulomb repulsion. Therefore, knowing the 
Pl value might only provide us with the instantaneous velocity of the fragments at 
this point, which could be obtained much more easily from conservation of energy. 
By contrast, it could be argued that the avoided crossing approaches the point of 
rupture as the fragment asymmetry increases. This suggests that the Pl value also 
becomes more sensitive to the pre-fission scenario in this case. 

Furthermore, in a series of papers |341 133 j it was noted that the point of rupture 
itself can be considered as an irregular point where the analyticity of the dependence 
of R(t) on time is broken. That is rather apparent from the mathematical viewpoint, 
as the analyticity would assume a unique trajectory in all the domain of its definition, 
< R < oo. Actually, the trajectory can be considered unique only after scission, 
for R sc < R < oo. But before scission, one can provide a rather arbitrary scenario, 
as we saw when discussing the one- and two-body mechanisms of friction. Moreover, 
the very position of the scission point itself may be varied within certain limits. We 
return to this question in a later section. 

In fact, the muon appears to be the only available tool for such studies. However, 
this simple picture is complicated by the fact that transitions to some of the higher- 
lying levels of the transient muonic molecule (e.g. 2pir and 2sa) result again in 
muon attachment to the heavy fragment. 

To obtain an order-of-magnitude estimate for the muon attachment probabilities, 
we utilize a simple formula derived by Demkov |36| and by Meyerhof |37j . Their 
model is based on the two lowest molecular levels (lscx and 2pcr) and utilizes first- 
order perturbation theory to calculate the transition probability from the lscr to the 
2p<7 level; within the two-level model, it is equal to the muon attachment probability 
to the light fission fragment: 



where Ih and II denote the binding energies of the muonic K-shell belonging to the 
heavy and light fission fragments, respectively, and v is the relative velocity of the 
fission fragments. For a fragment charge asymmetry £ = Zh/Zl = 55/39 = 1.41 
(corresponding to the peak of the mass distribution) one finds 1^=5.93 MeV and 
I_L=3.45 MeV. If we assume that the molecular transition occurs at a relative velocity 
of the fission fragments of v = 0.08c we obtain a muon attachment probability 
Pl = 0.042 for the light fission fragment; this value decreases to Pl = 0.015 if 
v = 0.06c. The choice of the parameters in the Demkov model and its extension are 
discussed in more detail in the next section. 




(3) 
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3 Physical premises 

3.1 Analytical properties of the trajectories 

The concept of determining the role of the avoided crossing for the transition prob- 
ability is a milestone in the quasiclassical theory of collisions accompanied with 
reorganization of the scattering systems, in the case where one or two of the sub- 
systems are described quasiclassically, and the other subsystems are described in 
terms of quantum mechanics. Most famous are collisions of the Landau-Zener type. 
Well-known are also the Demkov and Nikitin models. 

These models are based on the different character of the energy term's behavior 
in the vicinity of the scission point. In the two-level approximation, the usual 
expression for the terms can be written as follows: 

£i,2 = ^(ei + e 2 )±V / ( ei - e2 ) 2 + A2 • ( 4 ) 

The physical meaning of the parameters ei, e 2 in eq. (JIJ) is that they are consid- 
ered to be the non-interacting levels. The parameter A describes their interaction, 
and Ei, E2 are the true terms, with the interaction taken into account. 

The Landau-Zener model deals with intersecting terms e±, €2- The interaction 
A is then considered to be A = constant in the vicinity of the crossing. In the 
opposite extreme, the Demkov model considers the case of e\, e 2 parallel to each 
other, which is usually the case at large distances. At smaller distances, the terms 
start to diverge due to the interaction, which is supposed to be of exponential form: 

ei — 62 = e = constant, A = exp(-XR) . (5) 

The Nikitin model generalizes both the Landau-Zener and Demkov models with 
the following parameters: 

Cl - e 2 = e - Bcos9exp(-XR) , A = B sin 9 exp(-A-R) . (6) 

Specifically, the Demkov model is obtained with 9 = ir/2. The level crossing occurs 
in the complex plane at the point 

R c = (e/B)exp(i9) . (7) 

Many characteristic features of the muon distribution were consecutively stud- 
ied within the framework of the Demkov |38[ 140] and Nikitin models. The 
parameters of the models were related to realistic muonic wavefunctions. All the 
parameters, as well as the relative velocity, entering the expression for the transi- 
tion probability, have been taken at the point of the avoided crossing. Naturally, 
pre-fission dynamics was not involved in such an approach, as the pseudocrossing 
point is beyond the scission point as mentioned in the Introduction. Conclusions 
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about the pre-fission dynamics could be drawn from comparison with experimental 
data (e.g., [12]). Moreover, prospects of the study of other manifestations of the 
fission dynamics, such as muon shake-off due to the neck rupture, have been dis- 
cussed. These effects are not in the scope of the quasimolecular adiabatic picture 
considered herein. The results obtained were also applied for the interpretation of 
many experiments concerning muonic conversion of the prompt fission fragment 7 
rays, and muon capture by the fragments (e.g., jTH] and refs. cited therein). 



3.2 The complex trajectory method 

What was stated previously about the analytical properties of the trajectory function 
R(t) can be realized in terms of the complex trajectory method |34 | 1351 152] . 



1 3 



5 ~\ 



R. 



Rr 



Re R(t) 



Figure 1: The complex trajectory. The point R = corresponds to the united atom limit. The 
onset of population of the 2pa level is at the saddle point R s . The scission point R sc and the point 
of avoided crossing R c are the other singular points, where the integration contour is fixed. 

According to this method, the path R{t) corresponding to the trajectory of the 
relative motion of the fragments from R = to R — > oo can be deformed and shifted 
to the upper half-plane of complex R values, with the exception of irregular points. 

The reason for such a transformation is that in the upper half-plane the nonadia- 
batic matrix element vanishes. As a result, the main contribution comes from small 
segments nearest to the real axis of R, attached to the fixed irregular points. This 
transformation is shown in Figure^ The irregular points are, as was mentioned pre- 
viously, the pseudocrossing point i? c , which is a square-root type branching point, 
the scission point R sc , and the starting point of the motion which in our case is 
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reasonable to put on the top of the barrier. Thus, the contributions from the line 
segments 2 and 5 are expected to vanish for R — > oo |47j . 

We then note that the contributions from the segments 3 and 4 partly cancel 
one another, being of the opposite sign. Were the scission point a regular point of 
the trajectory, the cancellation would be complete, and this point would give no 
contribution to the transition amplitude. This would be the same as the point Rsc 
would be shifted to the upper half-plane, together with the segments 3 and 5. The 
two vertical cuts in Figure ^ fix the trajectory at the scission point, not allowing 
one to move it from the real axis of R. 

Relative contributions of the scission and pseudocrossing points turn out to be 
different, depending on the mass split of the fragments. The first contribution is 
several times smaller than the second one in the case of the most probable mass split. 
In the case of larger asymmetry, the contribution from scission approximately holds. 
In contrast, the contribution from the avoided crossing exponentially decreases. For 
this reason, the attachment probability becomes relatively more sensitive to the 
pre-fission scenario for more asymmetric fission, though its absolute value decreases. 
This circumstance was noted in refs. |46l 144] . where a higher sensitivity for larger 
asymmetry was predicted as a consequence of approaching the pseudocrossing point 
to the prefission area. 

We note that in actual calculations, the relative velocity V{R) is smoothly joined 
at the point R sc from the left and the right, with its derivative. Then the position 
of R sc is the only parameter which governs the contribution from the scission point. 
In this schematic picture, the pseudocrossing point appears to be completely inde- 
pendent of the prefission scenario, as the values of integrands along the segments 6 
and 7 are determined by analytical continuation of the Coulomb part of the trajec- 
tory at R sc < R < oo. This explains why in some papers the resulting transition 
probability was found not to depend on the prescission scenario: as a matter of fact, 
the position of R sc was not varied in those studies, but only the prefission time t sc 

Hang. 

3.3 The Born - Oppenheimer expansion 

Many of the first papers devoted to the calculation of the muon distribution were 
based on the Born - Oppenheimer expansion. We use the following series expansion 
0: 

V(r;R(t)) = Y J C n <&n{r]R{t))e i ^ r h- i P E '^' )dt ' . (8) 

n 

Here ^(r;R(t)) is the wavefunction, <3? n form a full set of quasimolecular wave- 
functions. Our basic product functions $ n (r; i?)e ip " r satisfy the time-dependent 
Schrodinger equation 

{i^-H)$ n {r-Ry™- iE '^ = 0, (9) 
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with modified eigenvalue E' n = E n + reflecting conservation of energy, and 

the two-centre Hamiltonian 

H = -^ + V 1 (\r-K 1 \) + V 2 (\r-R 2 \). (10) 

When the internuclear distance R — > oo each of the functions $ n correlates with 
a certain muonic state n of the corresponding fragment. 

The momentum-translation exponents (m.t.e.) are introduced in eq. (|SJ), where 
p n = fJ,v n , with n the muon mass, and v n a constant, equal to the velocity of 
the fragment in the asymptotic region. This condition is not essential, and can be 
easily abandoned [34 a . On the other hand, the orthogonality of the basis functions 
is broken in this description. 

Substituting the expansion (1) into the time-dependent Schrodinger equation, 
and neglecting small non-orthogonal terms, we get the following set of coupled equa- 
tions 

,IC ' ~J2F ik C k , (11) 



dR 



F ik = Mike^W-zlWVVWW, (12) 
M lk =< ^e l(Pir) l^ " ^^ V \^e t{pkr) > , (13) 

where the differential operators act only on the wavefunctions <3?j or & k , and V(R) 
stands for the relative velocity of the fragments. 

The equations are to be solved with the initial condition 

Ci = 5a for R = Ri , (14) 

denoting that the muon is in the ground state at the starting point 

Note, that for R — > oo, the modified basis regains its primordial orthogonality, 
as pi are the same for each state i on one fragment. Therefore, the amplitudes of 
decomposition over the modified basis acquire the usual physical meaning. These 
amplitudes must be used for the determination of physical probabilities, rather than 
the amplitudes in the traditional basis without the m.t.e. 

Incorporation of the m.t.e. thus modifies the non-adiabatic matrix element from 

m< - A (15) 

to 

M =dR + W(R)d-z ' (16) 
according to eq. (|13jl with the momenta p n = const. 
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Equation (|1U|) poses a two-dimensional eigenvalue problem. It has been solved 
numerically in ref.^lj by making use of the finite element method. On the other 
hand, it should be noted that using in eq. (jHJ) the asymptotic values of p n = const 
is only justified in the limit R — > oo. For small R, where the relative velocity of the 
fragments is several times smaller, the perturbation brought by the second terms in 
eqs. (j!3j) . (|16|) to the non-adiabatic transition operator is in principle excessively 
large. This difficulty was successfully overcome |34j . 

The problem (fTU|) has been solved in spheroidal coordinates 

^,77 = (n±r 2 )R (17) 
l<V<oo, -I<r/<1 . (18) 

It is noteworthy that recently, spheroidal coordinates were also applied in refs. 
|48| I49j , where the dynamics was studied of light charged particle emission in spon- 
taneous cold fission of 252 Cf. In these coordinates, the Laplacian becomes 

= w - v 2 ) W^ 2 ~ % + ^ (1 ~ v \ ] ■ (19) 



A finite-differences method on a rectangular grid (|18|) is used to obtain a first 
approximation to the solution of the problem (|lUj) . generally a rather poor one. 
The final solution of the problem is found by minimization of the Rayleigh-Ritz 
functional of the problem (|l(Jj) : 

where 

^) = |[W 2 -l)(||) 2 + (l-. 2 )(^) 2 ] , (21) 

and 

7? 3 

P(^V) = T 0> 2 -V 2 ) ■ (22) 

The minimization is achieved by the finite-element method. For realization of the 
method, all the area G of variables (|18|) is devided into N rectangular pieces g 3 . 
Within each element, a biquadratic function is defined as follows: 

(f>j (tp, rf) = a J + a{ ip + a\r\ + a^iprj + a\ ip 2 + a\ rf + a^i/j 2 !) + a^rf . (23) 

In the basis of functions <pj, j = 1,...N, the problem (|2*U)) is reduced to a generalized 
algebraic eigenvalue problem = EM&, with the "stiffness" matrix K, and M 
the matrix of "mass". The problem is then solved by the inverse iteration method 

Many results concerning the influence of the fission dynamics on the muon at- 
tachment probabilities were obtained by direct numerical integration of the time- 
dependent Schrodinger and Dirac equations on a 3-dimensional grid. This method 
and the results obtained are presented in the next sections. 
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4 Theoretical Developments 

Because the nuclear excitation energy in muon-induced fission exceeds the fission 
barrier height it is justified to treat the fission dynamics classically (no barrier 
tunneling). For simplicity, we describe the fission path by one collective coordinate 
R; the classical collective nuclear energy has the form 

£nuc = \b{R)R 2 + V &B (R) + E^R). (24) 

We utilize a coordinate dependent mass parameter B(R) [I] and an empirical double- 
humped fission potential Va s (R) which is smoothly joined with the Coulomb 
potential of the fission fragments at large R. The last term in Eq. (|24|) denotes the 
instantaneous muonic binding energy which depends on the fission coordinate; this 
term will be defined later. 




10 20 30 

R (fm) 



Figure 2: Phenomenological fission potential for 2 $lNp, for a mass asymmetry parameter 
Ah/Al = 1.40. The potential is constructed from four parabolic sections and is smoothly joined 
with the Coulomb potential of the fission fragments at large distance R. 

To account for the nuclear energy dissipation between the outer fission barrier 
and the scission point, we introduce a friction force which depends linearly on the 
velocity. In this case, the dissipation function D is a simple quadratic form in the 
velocity 

Enuc(t) = -2D = -fR 2 (t). (25) 

The adjustable friction parameter / determines the dissipated energy; it is the only 
unknown quantity in the theory. 
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For the dynamical description of the muonic wavefunction during prompt fission, 
the electromagnetic coupling between muon and nucleus (— e^f^A^) is dominant; the 
weak interaction is negligible. Because of the nonrelativistic motion of the fission 
fragments the electromagnetic interaction is dominated by the Coulomb interaction 

A°(r,t)= [ dV ^ 1 ^ . (26) 
J |r — r'| 

The muonic binding energy in the ground state of an actinide muonic atom amounts 
to 12 percent of the muonic rest mass; hence non-relativistic calculations, while qual- 
itatively correct, are limited in accuracy. Several theory groups have demonstrated 
the feasibility of such calculations [§1 El El which are based on the time-dependent 
Schrodinger equation 

h^V 2 - eA°(r,t)] VM) = ih^ i>{v,t). (27) 

In 1992, Oberacker et al. j!4| l4"| developed a numerical algorithm to solve the rela- 
tivistic Dirac problem on a three-dimensional Cartesian lattice. The time-dependent 
Dirac equation for the muonic spinor wave function in the Coulomb field of the fis- 
sioning nucleus has the form 

H D (t) ip(r,t) = ih^ if>(r,t), (28) 

where the Dirac Hamiltonian is given by 

H D (t) = -ihca ■ V + (3mc 2 - eA°(r, t). (29) 

Our main task is the solution of the Dirac equation for the muon in the presence of 
a time-dependent external Coulomb field A°(r,t) which is generated by the fission 
fragments in motion. Note the coupling between the fission dynamics, Eq. ([24 j). 
and the muon dynamics, Eq. (|28[). via the instantaneous muonic binding energy 

E ll (R(t)) = W(r,t)\H 13 (t)\il>(T,t)) (30) 

which depends on the fission coordinate; the presence of this term increases the 
effective fission barrier height. 

5 Lattice Representation: Basis-Spline Expansion 

For the numerical solution of the time-dependent Dirac equation (|28|l it is convenient 
to introduce dimensionless space and time coordinates 



x = r/A c A c = h/(m^c) = 1.87/m 
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r = t/r c t c = A c /c = 6.23 x 10~ M s (31) 

where A c denotes the reduced Compton wavelength of the muon and r c the reduced 
Compton time. For the lattice representation of the Dirac Hamiltonian and spinor 
wave functions we introduce a 3-dimensional rectangular box with a uniform lattice 
spacing Ax. The lattice points are labeled (x a ,yp, z^). 

Our numerical algorithm is the Basis-Spline collocation method jSZ]. Basis- 
Spline functions Bf^x) are piecewise-continuous polynomials of order (M — 1). 
These may be thought of as generalizations of the well-known "finite elements" 
which are B-Splines with M = 2. 




Figure 3: Set of Basis Splines of order M = 5 using periodic boundary conditions. 

There exists an extensive literature on Basis-Spline theory, developed by mathe- 
maticians |56j . In 1991, Umar et al. [57] applied the method to the numerical solution 
of problems in atomic and nuclear physics on a lattice; this paper discusses the B- 
Spline collocation method, periodic and fixed boundary conditions, the solution of 
various 1-D radial problems (Schrodinger, Poisson and Helmholtz equations), and 
the solution of 3-D Cartesian problems (Poisson equation). In a later paper in 1991, 
Umar et al. solved the nuclear static and time-dependent Hartree-Fock equations 
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on a 3-D lattice [SHI- During 1992-1995, Oberacker et al-OH] and Wells et al -Pl 
applied the B-Spline collocation method to the static and time-dependent Dirac 
equation which overcame the "fermion doubling problem" which one encounters 
with the finite-difference method. Very recently, we have used both the B-Spline 
collocation and Galerkin methods for the solution of the Hartree-Fock-Bogoliubov 
nuclear structure problem |6L)l I61j . 

To illustrate the B-Spline collocation method let us consider a wave function 
which depends on one space coordinate x; we represent the wave function on a finite 
spatial interval as a linear superposition of B-Spline functions 

TV 

^{x a )=Y J Bf 1 {x a )c\ (32) 
i=l 

In the Basis-Spline collocation method, local operators such as the EM potential A 
in Eq. Q29JI become diagonal matrices of their values at the grid points (collocation 
points), i.e. V{x) —* V a = V(x a ). The matrix representation of derivative operators 
is more involved [SH]. For example, the first-derivative operator of the Dirac equation 
has the following matrix representation on the lattice 

2V 

DP^B^B* , (33) 

i=i 

where B' ai = [dB^ {x) / dx\\ x=Xa . Furthermore, we use the shorthand notation B^ = 
B]^{xpt) for the B-spline function evaluated at the collocation point x», and the 
inverse of this matrix is denoted by B l @ = [B -1 ]^. Because of the presence of 
this inverse, the operator D@ will have a nonsparse matrix representation. In the 
present calculations we employ B-Splines of order M = 5. Eq. (|32|) can readily be 
generalized to three space dimensions; in this case the four Dirac spinor components 
ip(p\p = (1, • • 4) are expanded in terms of a product of Basis-Spline functions 

^(x a ,yp,z„t) = £ BV(x a )B?{v P )B?{zjj*{t), (34) 

i,3,k 

i.e. the lattice representation of the spinor wave function is a vector with N = 
A*N x *N y *N z complex components. Hence, it is impossible to store Hd in memory 
because this would require the storage of iV 2 complex double-precision numbers. We 
must therefore resort to iterative methods for the solution of the matrix equation 
which do not require the storage of Hj). 

We solve the time-dependent Dirac equation in two steps: first, we consider the 
static Coulomb problem at time t = 0, i.e. the muon bound to an actinide nucleus 



H D (t = 0)V> gs = E gs ip gs . 



(35) 
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This static problem is solved by an iterative procedure (damped relaxation method 
[3J02I]). The second part of our numerical procedure is the solution of the time- 
dependent Dirac equation (|2*%|) by a Taylor-expansion of the propagator. For an 
infinitesimal time step At we find 

w + At) = u(t + At, t)m « (i + y H#At)" 

We have thus reduced the original problem to a series of (matrix) x (vector) opera- 
tions which can be executed with high efficiency on vector or parallel supercomputers 
without explicitly storing the matrix in memory. 

6 Discussion of Numerical Results 

In the following we present results for prompt fission of 2 gjNp induced by the E2 : 
(3d — ► Is, 9.5MeV) muonic transition. All results reported here are for a 3-D Carte- 
sian lattice of size L x = L y = 67 fm and L z = 146 fm with N x * N y * N z = 25 * 25 * 53 
lattice points with a uniform lattice spacing Ax = 1.5A C = 2.8/m. Depending on 
the value of the friction coefficient, we utilize between 1, 200 — 1, 900 time steps with 
a step size At = 1.5r c = 9.3 x 10~ 24 s. The computer source code is written in 
Fortran 95. Typical production runs use 25 MB of memory and require 5.3 hours 
of CPU time on an Intel Pentium-4 (1.7GHz) GNU/LINUX workstation. The same 
job takes 3.2 CPU hours on an IBM-SP2 supercomputer at NERSC (serial run on 
one processor). 




Figure 4: Prompt muon-induced fission of ^Np for a fission fragment mass asymmetry £ = 
Ah/Al = 1.10 at E* — 9.5 MeV. Shown is the Coulomb interaction energy between the muon and 
the fissioning nucleus (left) and the corresponding muon position probability density (right) at time 
t = during fission. Zero friction (/ = 0) has been assumed. 
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Figure 5: Same as Fig. 4, except at time step f = 6.5x 10 21 s. 




Figure 6: Same as Fig. 4, except at time step f = 8.4x 10 21 s. 

Fig. 0] shows the Coulomb interaction energy between the muon and the neptu- 
nium nucleus (left) and the corresponding muon position probability density (right) 
at time t = during fission. At this moment, the actinide nucleus is situated in its 
deformed ground state minimum with quadrupole deformation (3% = 0.27. Figures 
!5l6l7l show the time-development of the Coulomb interaction and the corresponding 
muon position probability density during fission. The calculation has been done for 
an asymptotic fragment mass asymmetry £ = Ah/Al = 1.10. In the last two time 
steps the Coulomb potential wells of the two separated fission fragments are clearly 
visible; the deeper well on the left is generated by the heavier fission fragment. We 
observe a clear spatial correlation between the positions of the potential minima 
and the muon probability density maxima. As expected, the muon sticks predomi- 
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nantly to the heavy fragment (large bump on the left in Fig(7J), but since the mass 
asymmetry in this case is very small (£ = Ah/Al = 1.10) the muon attachment to 
the light fission fragment (bump on the right) is relatively large. By integrating the 
muon probability densities associated with the heavy and light fission fragments at 
large internuclear distances, we can infer the muon attachment probabilities to the 
heavy and light fragments, Ph and Pl, respectively. From Fig(7|we find Pl = 0.20. 

One might ask whether the muon will always remain bound during fission; what 
is the probability for ionization? To investigate this question we have plotted the 
muon position probability density on a logarithmic scale. 




10 20 30 40 50 



Figure 8: Contour plot of the logarithm of the muon probability density at t = 1.1 x 10 20 s 
shows no evidence of muon ionization. 

In coordinate space, any appreciable muon ionization would show up as a "prob- 
ability cloud" that is separating from the fission fragments and moving towards the 
boundaries of the lattice. Fig. |H] shows no evidence for such an event in our numer- 
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ical calculations. Hence, we conclude that the probability for muon ionization P[ on 
is substantially smaller than the muon attachment probability to the light fission 
fragment which is always clearly visible in our logarithmic plots, even at large mass 
asymmetry. From this we estimate that Pi on < 10~ 4 . 

As in all lattice calculations, we need to demonstrate convergence of our results 
in terms of the lattice size and lattice spacing. Fig. |§] shows the asymptotic muon 
attachment probability Pl as a function of the total number of lattice points. In all of 
our production runs we use a lattice with a total of N x *N y *N z = 25*25*53 = 33, 125 
lattice points. It is apparent that the muon attachment probabilities (and other 
related observables) have indeed converged. 
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Figure 9: Demonstration of numerical convergence of muon attachment probability for 2 gjNp. 
Calculations use a mass asymmetry of 1.10 and zero friction. 

Fig-E3sh° ws that Pi depends strongly on the fission fragment mass asymmetry. 
This is easily understood: for equal fragments we must have Pl = 0.5, and for large 
mass asymmetry it is energetically favorable for the muon to be bound to the heavy 
fragment, hence Pl will be small. 

Fig. llll shows the nuclear energy dissipation (in form of neutron- and 7-emission) 
as a function of time; in our model, friction is confined to the region between the 
outer fission barrier and the scission point; for friction parameters / = 200 and / = 
500, we obtain total dissipated energies E^ sa = 14.2MeV and 22.0MeV, respectively. 

In Fig. we examine the dependence of Pl on the dissipated nuclear energy, 
i?diss> during fission. In our model, friction takes place between the outer fission 
barrier and the scission point. When the dissipated energy is computed from equa- 
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Figure 10: Muon attachment to light fission fragment vs. fission fragment mass asymmetry, 
upper curve: 2 ff U, lower curve: 2 f|Np. The calculations assume zero friction. 
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Figure 11: Energy dissipated during fission, for friction parameters / = 0, / = 200 and / = 500. 



tion (|25|) we find an almost linear dependence of the muon attachment probability 
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on -Edissj unfortunately, this dependence is rather weak. 




5 10 15 20 25 



E diss ( MeV ) 

Figure 12: Muon attachment probability to the light fission fragment as function of nuclear 
energy dissipation for 2 ||Np. Results are shown for fragment mass asymmetries £ = 1.05 (upper 
curve), 1.10, 1.15, and 1.20 (lower curve). 

We would like to point out that the theoretical values for Pj, obtained in this 
work are smaller than those reported in our early calculations |14|l4"j. There are two 
reasons for this: (a) the size of the lattice and (b) the lattice representation of the 
first derivative operator in the Dirac equation. Because of constraints in the amount 
of computer time available to us we utilized a smaller cubic lattice in our prior 
calculations [3] with N x * N y * N z = 29 3 lattice points. In recent years, we were able 
to increase the size of the lattice substantially, in particular in fission {z-) direction 
(see Fig. Regarding the lattice representation of the first derivative operator, 
Eq. (|33|) . in the Dirac equation: in ref. |141 H] we utilized a combination of forward 
and backward derivatives for the upper and lower spinor wave function components; 
after extensive testing of Coulomb potential model problems with known analytical 
solutions we have found that the symmetric derivative operator provides a more 
faithful lattice representation. The results reported here and in ref. |15j have been 
obtained utilizing the symmetric derivative prescription. 

7 Comparison of Theory with Experiment 

Prompt muon-induced fission was first observed experimentally by Diaz et al. |62j . 
More recent experiments by Ahmad et al. ^5] yield a total fission probability per 
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muon stop Pf = 0.068 for 238 U and a ratio Pf (prompt) /Pf (delayed) = 0.089. 

There are only few experimental data on muon attachment available for compar- 
ison with our theory. Schroder et al. 181 measured for the first time mean lifetimes 
of muons bound to fission fragments of several actinide nuclei. The muon decays 
from the K-shell of the muonic atom through various weak interaction processes at 
a characteristic rate A = Ao + A c , where Ao = (2.2 x 10~ 6 s) -1 is the free leptonic 
decay rate for the decay process pT — > e~ + v e + and A c denotes the nuclear 
capture rate; A c depends upon the charge and mass of the fission fragment. From 
the observed lifetime = 1.30 x 10~ 7 s Schroder et al. estimated an upper limit 
for the muon attachment probability Pl < 0.1. It must be emphasized that this 
number represents an integral over the whole fission mass distribution and, hence, 
cannot be directly compared to the numbers given in Fig. 1121 

The most complete experiments have been carried out by Risse et al. us- 
ing the muon beam of the ttE3 channel at the Paul Scherrer Institute (PSI) in 
Switzerland. For this purpose, a fission chamber has been inserted into the elec- 
tron spectrometer SINDRUM I. The incident muons are detected by a scintillation 
counter. An event is defined by a {fi>~ , fif2&~) coincidence where the fission frag- 
ments are observed in prompt and the muon decay electrons in delayed coincidence 
with respect to the incident muon. The magnetic field of the electron spectrometer 
allows for a reconstruction of the electron trajectories. Thus, it is possible to deter- 
mine whether the muon decay electrons originate from the heavy or the light fission 
fragment. 

For several mass bins of the light fission fragment, muon attachment probabili- 
ties Pl have been measured; the experimental data are given in Tabled R should 
be emphasized that the mass bins are relatively broad. Because the theoretical 
values for Pl depend strongly on the mass asymmetry it is not justified to assume 
that Pl remains constant within each experimental mass bin. Instead, to allow for 
a comparison between theory and experiment, we have to multiply the theoretical 
Pl values in Fig. El with a weighting factor that accounts for the measured relative 
mass distribution of the prompt fission events within this mass bin. We subse- 
quently integrate the results over the sizes of the experimental mass bins. Due to 
the relatively low excitation energy in muon-induced fission, the fission mass dis- 
tribution exhibits a maximum at £ = Ajj/Al = 1.4 and falls off rather steeply for 
values larger or smaller than the maximum. This means that the large values of 
Pl ~ 0.5 at or near fission fragment symmetry £ = 1.0 will be strongly suppressed. 
The resulting theoretical values for Pl are given in the last column of Table H R 
is apparent that our theory agrees rather well with experiment. Because of the 
size of the error bars in the experiment and because of the weak dependence of the 
theoretical values of Pl on the dissipated energy, it is not possible to extract very 
precise information about the amount of energy dissipation. 

From a comparison of our theoretical result for the mass bin Al = 118.5 — * 111.5 
with the measured data we extract a dissipated energy of order 10 MeV for 237 Np 
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Table 1: Muon-attachment probabilities to the light fission fragment, Pl, for 237 Np(/i , /). Exp. 
data are taken from ref.|25|. 



mass bin Al 


mass asymmetry 


P L (exp) 


P L (theo) 


118.5 -> 111.5 


1.000 -> 1.126 


(25.5 ±8.5) x 10~ 2 


26.0 x 10- 2 ,E diBB = OMeV 
22.3 x 10 _2 ,-B diss = 22MeV 


111.5 -> 104.5 


1.126 -> 1.268 


(9.7 ±2.6) x 10" 2 


6.62 x W- 2 ,E diBB = OMeV 
3.51 x 10" 2 ,£ diss = 22MeV 



while the second mass bin Al = 111.5 — * 104.5 is more compatible with zero 
dissipation energy. We like to point out that the value E^ iss = 10 MeV agrees 
with results from other low-energy fission measurements that are based on the odd- 
even effect in the charge yields of fission fragments |65| . In addition to 237 Np we 
have also studied muon-induced fission of 238 U; the results for muon attachment are 
very similar (see Fig. llOj) . 

8 Outlook: calculation of dissipation in fission 

By making use of the commutators [^,H] , one can express the non- 

adiabatic matrix elements (|16|) in terms of the mean-field derivatives and the energy 
differences of the initial and final energies: 

„ 1 ( d , v x - v 2 d 



^ = -s^\M + wmBi) m ■ <37> 

V(r) in eq. (|37jl is the mean field potential. In the two-center model it can be 
written as follows: 



Vx(\f-Rx\) forz<0 
Hr) J (38) 
V 2 (\f-R 2 \) forz>0 

As the derivative is essentially non-zero at the nuclear surface, this is in the 
spirit of the well-known semi-classical wall-and-window model, where dissipation 
arises from the interaction with the walls of a fissile nucleus. 

We note that the change of the mean potential ()38|) during separation of the 
fragments is due to two reasons: One is the change of the potential due to its radial 
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dependence through Vi(\r — Ri\). The corresponding contribution to the matrix 
element (|57|) cancels |E] due to the equality 

as*!-. <*> 

which can be obtained by direct calculation. 

The other contribution is due to "breathing" of each of the nascent fragments, 
which arises from the change of its nuclear radius Ro: 

% M) _««£z^* . (4 „) 



dRi Vl u dR dR 

In the form (|4()|) . the perturbation Hamiltonian was applied for calculating the 
energy dissipated during saddle to scission descent j^j within the two-center har- 
monic oscillator model. The calculated amount turns out to be rather small (about 
1 MeV) which gives evidence for a rather undamped motion for this kind of friction. 
The result is in agreement with what was said in the Introduction about this mech- 
anism of nuclear friction. It also agrees with the conclusion drawn from comparison 
of the calculated muon attachment probabilities with experiment. 



9 Conclusions 

In this paper, we study the process of prompt muon-induced fission of actinide nuclei. 
The nuclei are excited by radiationless transitions (inverse internal conversion). For 
example, in 2 || U the El : (2p — ► Is) and the E2 : (3d — ► Is) muonic transitions 
result in excitation of the nuclear giant dipole and giant quadrupole resonances, 
respectively, which act as doorway states for fission. It is very probable that the 
(3d — * Is, 9.5MeV) transition in the muonic atom will be dominant for muon- 
induced fission, because the transition energy is very close to the peak of the T = 
giant quadrupole resonance of the nucleus (Eqqr = 9.9MeV). By contrast, there 
is a mismatch between the (2p — > Is, 6.6MeV) muonic transition energy and the 
center of the giant dipole resonance which is located at Eqdr = 12.8MeV. Because 
the muon lifetime is long compared to the timescale of prompt nuclear fission, the 
motion of the muon in the Coulomb field of the fissioning nucleus may be utilized 
to learn about the dynamics of fission. 

We have studied the dynamics of a muon bound to a fissioning actinide nucleus 
in two different ways: first, using the non-relativistic time-dependent Schrodinger 
equation, the Born-Oppenheimer expansion method was utilized. This method relies 
on a set of quasimolecular wavefunctions, and the fission dynamics is described 
classically. 

In the second method, we have solved the relativistic time-dependent Dirac equa- 
tion on a 3-D Cartesian lattice using Basis-Spline expansion techniques. (The 3-D 
lattice makes it possible to investigate nuclear fission shapes that are not axially 
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symmetric. However, in the present study we always assume axial nuclear sym- 
metry, for simplicity). The principal advantage of the coordinate lattice method is 
that it does not depend on the validity of an "atomic" or "molecular" wavefunc- 
tion basis set. Furthermore, the method accounts not only for transitions between 
bound muonic states, but is also describes the possible ionization of the muon dur- 
ing fission: in coordinate space, any appreciable muon ionization would show up 
as a "probability cloud" that is separating from the fission fragments and moving 
towards the boundaries of the lattice. Fig. |H] shows no evidence for such an event in 
our numerical calculations. 

Our time-dependent Dirac equation calculations predict a strong mass asymme- 
try dependence of the muon attachment probability Pl to the light fission fragment; 
this feature is in agreement with experimental data. The theory also predicts a (rel- 
atively weak) dependence of Pl on the dissipated energy. By comparing our theo- 
retical results to the experimental data of ref. we extract a dissipated energy 
of order — 10 MeV for 237 Np (see Table 1). Using the dissipation function defined 
in Eq. (|25jl. the 10 MeV value corresponds to a fission time delay from saddle to 
scission of order 2 x 10~ 21 s. 
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